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We classify all spherically symmetric spacetimes admitting a kinematic self-similar vector 
of the second, zeroth or infinite kind. We assume that the perfect fluid obeys either a 
polytropic equation of state or an equation of state of the form p = K/i, where p and fi are 
the pressure and the energy density, respectively, and AT is a constant. We study the cases 
in which the kinematic self-similar vector is not only "tilted" but also parallel or orthogonal 
to the fluid flow. We find that, in contrast to Newtonian gravity, the polytropic perfect- fluid 
solutions compatible with the kinematic self-similarity are the Friedmann-Robertson- Walker 
solution and general static solutions. We find three new exact solutions which we call the 
^vq ' dynamical solutions (A) and (B) and /l-cylinder solution, respectively. 



■ §1. Introduction 

(N ' 

There is no characteristic scale in Newtonian gravity or general relativity. A 
set of field equations is invariant under a scale transformation if we assume appro- 
i , priate matter fields. This implies the existence of scale- invariant solutions to the 

field equations. Such solutions are called self-similar solutions. Among them, the 
spherically symmetric self-similar system has been widely researched in the context 
of both Newtonian gravity and general relativity. Although self-similar solutions are 
only special solutions of the field equations, it has often been supposed that they 
play an important role in situations where gravity is an essential ingredient in a 
spherically symmetric system. In particular, a self- similarity hypothesis has been 
proposed, which states that solutions in a variety of astrophysical and cosmological 
situations may naturally evolve to a self-similar form even if they are initially more 
complicatedliF . 

Self-similar solutions in Newtonian gravity have been studied in an effort to 
obtain realistic solutions of gravitational collapse leading to star formation 
For an isothermal gas cloud, Larson and Penston independently found a self-similar 
solution, which called the Larson-Penston solution, describing a gravitationally 
collapsing sphere BP' a>. Thereafter, Hunter found a new series of self-similar solutions, 
and noted that a set of such solutions is infinite and discrete &. Recent numerical 
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simulations and mode analyses showed that the Larson-Penston so^u^iqn gi^ 
best description for the central part of a generic collapsing gas sphereli 
For a polytropic equation of state, Yahil found the polytropic counterpart of the 
Larson-Penston solution describing a gravitationally collapsing sphere i3 (see also 
Refs. ||),[l3|)). We hereafter refer to these solutions as the polytropic Larson-Penston 
solutions. 

In general rel ativ ity, self-similarity is defined by the existence of a homothetic 
Killing vector fielded. Such self-similarity is called the first kind (or homothety). Ori 
and Piran discovered the general relativistic counterpart of the Larson-Penston self- 
similar solution together with Hunter's family of self-similar solutions for the perfect 
fluid obeying an equation of state p = K/s (0_< K < 0.036), where p and fi are the 
pressure and the energy density, respectivelyEj)'E3)'HZP. They observed that a naked 
singularity forms in this solution for < K < 0.0105. Harada and Maeda found that 
generic non-self-similar spherical collapse converges to the general relati vis tic Larson- 
Penston solution in an approach to a singularity for < K < 0.036t5''Ej). Since 
a naked singularity forms for < K < 0.0105, this implies the violation of cosmic 
censorship in the spherically symmetric case (see also Refs. |20|),|2T|)). This represents 
the strongest known counterexample against cosmic censorship. It also provides 
strong evidence for the self-similarity hypothesis in general relativistic gravitational 
collapse. The question then naturally arises, whether collapsing self-similar solutions 
with a polytropic equation of state exist in general relativity. If such solutions do 
exist, they may play an important role in the final stage of generic collapse, as in 
the p = Kfj, case. 

In Newtonian gravity, self-similarity for the polytropic case has a different form 
of the dimensionless variable from that in the isothermal case, since sound speed is 
not constant in the former case. The self-similarity coordinate is given by t 2 ~ 7 /r 
for the polytropic case, and by t/r for the isothermal case. The scaling functions 
of physical quantities for the former case are also different from those for the latter 
caseH3). In general relativity, there exists a natural generalization of homothety 
called kinematic self- similarity, which is defined by the existence of a kinematic self- 
similar vector fieldil) (see also the earlier related works by Tomita©). Kinematic 
self-similarity is characterized by an index and classified into three kinds: the second, 
zeroth and infinite kinds. 

One can show that an equation of state of the form p = K[i is the only barotropic 
one compatible with self-similarity of the first kind0\ Self-similar perfect-fluid so- 
lutions of the first kind with an equation of state of this form have been classified 
for the dust case (K = 0) by CarrcfP and for the case < K < 1 by Carr and 
ColeylH^ (see also Ref. ^6|)). Special cases in which a homothetic Killing vector 
is not "til ted" (i.e., either parallel or orthogonal to the fluid flow) have also been 
studiedc2MiiP. Kinematic self-similar perfect-fluid solutions have been explored by 
several authors ©'0' ©'11). Benoit and Coley have studied spherically symmet- 
ric spacetimes which admit a kinematic self-similar vector of the second and zeroth 
kinds©. Sintes, Benoit and Coley have considered spacetimes which admit a kine- 
matic self-similar vector of the infinite kindlSIP . In these works, the equation of state 
has not been specified. We have previously investigated spherically symmetric space- 
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times which contain a perfect fluid obeying a relativistic polytropic equation of state 
and admit a kinematic self=similar vector of the second kind in which the kinematic 
self-similar vector is tiltedS. There, we assumed two kinds of polytropic equation 
of state in general relativity and showed that such spacetimes must be vacuum in 
both cases. Although a spherically symmetric spacetime which contains a relativistic 
polytropic perfect fluid is incompatible with kinematic self-similarity of the second 
kind, it could be compatible with other kinds of kinematic self-similarities (i.e., the 
zeroth or infinite kind), or with the case in which a kinematic self-similar vector is 
parallel or orthogonal to the fluid flow. 

In this paper, we extend our previous work in several important ways. We study 
spacetimes which contain a perfect fluid obeying either a polytropic equation of state 
or an equation of state p = Kfi, and which admit a kinematic self-similar vector field 
of the second, zeroth or infinite kind. We assume two kinds of relativistic polytropic 
equations of state and study the case in which a kinematic self-similar vector is not 
only tilted but also parallel or orthogonal to the fluid flow. 

The organization of this paper is the following. In §|2| basic equations in a spher- 
ically symmetric spacetime are presented and kinematic self-similarity are briefly 
reviewed. We treat the cases in which a kinematic self-similar vector is tilted, par- 
allel and orthogonal to the fluid flow in §|3[ §||] and §||, respectively. §^ is devoted 
to summary and discussions. Possible equations of state which are compatible with 
self-similarity are discussed in appendix A. We adopt units such that c = 1. 

§2. Spherically Symmetric Spacetime and Kinematic Self-Similarity 

The line element in a spherically symmetric spacetime is given by 

where dfi 2 = d6 2 + sin 2 Qdtp 2 . We consider a perfect fluid as a matter field 



T^u = P(t, r) 9lxv + [n{t, r) + p(t, r)]U^U v , (2-2) 
U„ = (-e # , 0,0,0), (2-3) 

where is the four- velocity of the fluid element. We have adopted the comoving 
coordinates. Then the Einstein equations and the equations of motion for the perfect 
fluid are reduced to the following simple form: 



Pr 
fi+p' 

jit 2Rt 
H + P R 

iTT/lRrR 2 , 

-4irpR t R 2 , 



(2-4) 

(2-5) 

(2-6) 
(2-7) 
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-Rtr+$rRt + &tRr, 

1 



m 



2G 



Ril + e-^R^-e-^Rr 2 ), 



(2-8) 
(2-9) 



where subscripts t and r mean the derivative with respect to t and r, respectively 
and m(t, r) is called the Misner-Sharp mass. We also write the auxiliary equations: 



7^ 



Rf \ 2 Rf _ 
— +2—& t 
R R 



+ e 



20-2^ 



, Rrr 

'IT 



2'1> 



SirGfie 

£ I 2&-2<P 
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ir)' + 4* 



-8vrGpe 



2<f 



(2-10) 



(2-11) 



-20 



Rt 



Rt*t Rt$ 
R 



R 



(2-12) 



Five of the above nine equations are independent. 

In this paper, we assume the following two kinds of polytropic equations of state. 
One is 



p = K^T 1 , 

where K and 7 are constants, and the other isE 



(2-13) 



P 



p = m^n + 



V 



7-1' 



(2-14) 



where the constant m& and n(t, r) correspond to the mean baryon mass and the 
baryon number density, respectively. We call the former equation of state (I) (EOS 
(I)) and the latter equation of state (II) (EOS (II)). Here we assume that K 7^ and 
7 ^ 0, 1. We also treat an equation of state 



p = Kfi. (2-15) 

We call this one equation of state (III) (EOS (III)) and assume that — 1 < K < 1 
and K / 0. 

We note the properties of EOS (I) and (II). For 7 < 0, the fluid suffers from 
thermodynamical instability. For < 7 < 1, both EOS (I) and (II) are approximated 
by a dust fluid in high-density regime, since p/p = Kp 1-1 — > for \i — > 00 for EOS 
(I) and 
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~ — ^ Kn<-^ ~* for ^ oo, (2-16) 



7-1 



for EOS (II). For 1< 7, EOS (II) is approximated by EOS (III) with K = 7 - 1 in 
high-density regime since 

Kri^ Kri^ p 

fi = m&n H > = for n — > 00. (2-17) 

7 — 1 7 — 1 7 — 1 

In the case of 2 < 7 for EOS (II) and 1 < 7 for EOS (I), the dominant energy 
condition can be violated in high-density regime, which will be unphysical. 

In this paper, we introduce a vector field which is called a kinematic self-similar 
vector field and treat three cases in which it is parallel, orthogonal to the fluid flow 
and tilted, i.e., neither of them. In a spherically symmetric spacetime, a vector field 
7] is written in general as 

where h\(t,r) and h,2(t,r) are arbitrary functions. /12 = when 77 is parallel to the 
fluid flow, while h\ = 0, when r\ is orthogonal to the fluid flow. When rj is tilted, 
both hi and /12 are non-zero. 

A kinematic self-similarity vector ^ satisfies the condition 



£ 5 /V = 25V> (2' 19 ) 
= aUp, (2-20) 

where h^ v = g^ v + U^U U is the projection tensor, denotes the Lie differentiation 
along £ and a and 5 are constants E2rEiP. The similarity transformation is charac- 
terized by the scale-independent ratio, ot/6, which is referred to as the similarity 
index. 

In the case of 5 7^ 0, 5 can be set to unity and the kinematic self-similar vector 
£ can be written as 



«^ = ( < rt +%+'*- (2 ' 21) 

if it is tilted. In the case of a = 1, which is corresponding to self-similarity of the 
first kind {(3 can be set to zero), it follows that £ is a homothetic vector and the 
self-similar variable £ is given by r/t . In the case of a = 0, which is corresponding 
to self-similarity of the zeroth kind (fi can be rescaled to unity), the self-similar 
variable is given by 



£ = re 1 . 



(2-22) 
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In the case of a 7^ 0,1, which is corresponding to self-similarity of the second kind 
{(5 can be set to zero), the self-similar variable is given by 



f = , \ 1 / • (2-23) 

(at) 1 /" v ' 

In the case of 5 ^ 0, self-similarity implies that the metric functions can be written 
as 



R = rS(0, # = # = (2-24) 

In the special case of <5 = and the self-similarity is referred to as the infinite 

kind (a = 1 is possible). The kinematic self-similar vector £ can be written as 

e^- = 4+r^. (2-25) 
ox^ at or 

The self-similar variable is given by 



£ = - f (2-26) 

In the case of 5 = 0, self-similarity implies that the metric functions can be written 

as 



R = S(0, e* = e # «\ e* = e*«tyr. (2-27) 

It is noted that £ is a Killing vector in the case of 5 = and a = 0. If the kinematic 
self-similar vector £ is parallel to the fluid flow, it is written as 

d d 

= f(t)lZ, (2-28) 
s dx» w dt K J 

where f(t) is an arbitrary function and the self-similar variable is then r. If the 
kinematic self-similar vector £ is orthogonal to the fluid flow, it is written as 

(2 ' 29) 

where g(r) is an arbitrary function and the self-similar variable is then t. 

Mathematical and physical properties of spacetimes which admit a kinematic 
self-similar vector of the "finite" kind, i.e., tha second and zeroth kinds, and contain 
a perfect fluid have been discussed by Coleyca). It is noted that they have assumed 
"physical" self-similarity, i.e., 
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C^fi = a/i, C^p = bp, (2-30) 

where a and b are constants, in addition to kinematic self-similarity as represented 
by equations (2-19) and (2-20). The spacetimes which admit a kinematic self-similar 
vector of the infinite kind and contain a perfect fluid have been studied by Sintes, 
Benoit and Coley©. Benoit and Coley have investigated spherically symmetric 
spacetimes which admit a tilted kinematic self-similar vector field of the seconHor 
zeroth kind and contain a perfect fluid without specifying the equation of statelet*. 



§3. Tilted case 



3.1. Self- similarity of the the second kind 

In the case of self-similarity of the second kind, the Einstein equations imply 
that the quantities m, fj, and p must be of the form 



2Gm = r 



M 1 (C) + -M 2 (0 



8nGfi = 



8nGp 



(3-1) 
(3-2) 
(3-3) 



where £ = r/(at) l l a . A set of ordinary differential equations is obtained when one 
demands that the Einstein equations and the equations of motion for the matter 
field are satisfied for the 0[(r/t)°] and 0[(r/t) 2 ] terms separately. The equations for 
a perfect fluid (2-4)-(2Tl) reduce to the following: 



Mi + M[ 


= W!S 2 (S + S'), 


(3-4) 


3M 2 + M' 2 


= W 2 S 2 (S + S'), 


(3-5) 


M[ 


= -PiS 2 S', 


(3-6) 


2aM 2 + M' 2 


= -P 2 S 2 S', 


(3-7) 


Mi 


= S[l-e- 2 *(S + S')% 


(3-8) 


a 2 M 2 


= SS' 2 e- 2<p , 


(3-9) 


(Pi + Wx)®' 


= 2Px-P{, 


(3-10) 


(P 2 + W 2 )& 


- P' 


(3-11) 


W{S 


= -(P 1 + W 1 )(V'S + 2S / ), 


(3-12) 


(2aW 2 + W 2 )S 


= -(P 2 + W 2 )(1r'S + 2S'), 


(3-13) 


S" + S' 


= S'<P' + {S + s'p', 


(3-14) 


S'(S' + 2W'S) 


= a 2 W 2 S 2 e 2 *, 


(3-15) 
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2S(S" + 2S') - 2&'S(S + S') 
2S{S" + aS' - <P'S') + S' 2 
(5 + S')(S + S' + 2<P'S) 



i/2 



-s 

-a 2 P 2 S z e 
{l + P 1 S 2 )e 2 * 



S 2 + e 2 ^(l-W 1 S 2 ), 



(3-16) 
(3-17) 
(3-18) 



where the prime denotes the derivative with respect to ln£. 

In a vacuum case, while the Minkowski spacetime is compatible with all a, the 
Schwarzschild spacetime is compatible only with a = 3/2, since equations (3-5) and 
(3-7) are degenerated so that M 2 ^ is possible only in this case. The Schwarzschild 
spacetime in the Lemaitre's choice of coordinates is written as 



/ 



-dt 2 + r 2 J 3 



dp 2 



(P-t) 



2/3 



[3/ 


4/3 al 


- 2 (p-t)_ 


dn 2 



(3-19) 



where r g is a constant and the Schwarzschild radius corresponds to r g = (3/2) (p — 
i)EP. Changing the radial coordinate as p = r 3/2 , the metric in the form of (2-24) 
for a = 3/2 can be obtained: 



ds 2 = -dt 2 + rr 2/3 + r 2 



'(1-t/r 3 / 2 ) 



4/3 



df2 2 . (3-20) 



3.1.1. EOS (I) and (II) 

Subtracting equation (3-18) from equation (3-16) and eliminating S" with using 
equation (3-14), we obtain 



2<f>' = (P 1 + W{]e m . 
Then equations (3-10) and (3-11) result in 



(3-21) 



e^iPi + Wx) 2 = AP l -2P[, (3-22) 
e 2 *(P x + W 1 )(P 2 + W 2 ) = -2P' 2 . (3-23) 

If a perfect fluid obeys EOS (I) for K 7^ and 7/ 0,1, we find from equations 
(3-2) and (3-3) that 

a = 7 , P 1 = W 2 = 0, P 2 = (8 ^_i 72 r 27 ^i 7 , [case (A)] (3-24) 

or 

a=i P 2 = W 1 = 0, Pl = -*i—i*wq. [case(B)] (3-25) 

If a perfect fluid obeys EOS (II) for K 7^ and 7 / 0,1, we find from equations 
(3-2) and (3-3) that 
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a = 7 , P = 0, P 2 = m 7 (87rG)7 -i 72 r 27 ^? = (7 ~ 1)W 2 , [case(C)] 

(3-26) 

or 

« = i P = 0, Pi = m 7 (87r g )7 -i 727 e 2 ^ 2 7 = (7 - mi- [case (D)] (3-27) 

Pi = directly implies Wi = from equation (3-22), while Pj = implies 
(Pi + Wi)W 2 = from (3-23), which results in W 2 = for cases (B) and (D). 
Therefore the spacetime must be vacuum for all possible cases. 

3.1.2. EOS (III) 

If a perfect fluid obeys EOS (III), we find from equations (3-2) and (3-3) that 



P 1 = KWi, P 2 = KW 2 . [case (E)] (3-28) 

K ^ — 1 can be concluded, since K = — 1 implies Pi = W\ = from equations (3-10) 
and (3-12) and also P 2 = W 2 = from equations (3-11) and (3-13), which means a 
vacuum spacetime. For K ^ — 1, W\W 2 = can be proved. Assuming that W\ ^ 
and W 2 / 0, we find from equations (3-12) and (3-13) that W{/Wi = 2a + W£/W 2 , 
while equations (3-10) and (3-11) give that 2 - W[/Wi = -W 2 /W 2 . These two 
equations contradict the assumption of a ^ 1. Hence W\W 2 = is concluded. 

In the case of W\ = P\ = and IV2 / 0, M\ = can be obtained from equations 
(3-4) and (3-6), and then equation (3-8) gives (S + S") 2 = e 2 * '. Equation (3-21) gives 
<P' = 0, and then equation (3-11) requires P 2 to be constant, which implies that W 2 
is also constant. We set = cq and W 2 = wq, where cq and wq are constants. 
Equations (3-5), (3-7) and (3-9) give the evolution equation for S: 



3-2a{S'\ 2 M _ S' 

- (1 + K)w — - w = 0. 



2 n 2 



a*c; 



(3-29) 



The solution is S = so£ g , where sq and q are constants, q = — 2a/[3(l + K)] can be 
obtained from equation (3-13). The resulting solution is 



e = c , 



e = s 



1 



2a 



3(1 + K) 



2a 



3(1+K) 



S = SqC^ 1 ^, 

Mi = 0, M 2 
P 1 = W 1 = 0, 
P 2 = KW 2 = 



4ag 



9c 2 (l + ^) 2 
4if 



2a 



3c 2 (l + K) 2 ' 



(3-30) 
(3-31) 

(3-32) 
(3-33) 
(3-34) 
(3-35) 
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This is the flat Friedmann-Robertson- Walker (FRW) solution. Equation (3-31) im- 
plies a / 3(1 + K)/2. With the coordinate transformations if = c$t and r' = 
{cq/ , Q)~ 2 /[ 3 ( 1 +- ftr )] s^r 1 " 3 a+*o } the solution in a more usual form 

ds 2 = -dt' 2 + t'^^(dr' 2 + r' 2 dQ 2 ), (3-36) 

2a " = WTW^K- (3 ' 37) 

SxGp = 8*GJf„ = H1 * K K)V , (3-38) 

is obtained. 

In the case of W 2 = P 2 = and W\ / 0, equations (3-12) and (3T5) give 

?f£ + T^=0. (3-39) 



Eliminating 5" — S'<P' and "P 7 from equation (3-17) with using equations (3-14) and 
(3-12), respectively, we obtain 



S' / S"\ /.S" 1 W{ 



(a-1)-- 1 + - 2- + 



S V 5/ V S l + KWx 



S' 



2 



+ (-) =0. (3-40) 



If 3S'/S + [2/(1 + K))W[/W l = 0, equation (3-40) results in (2a - 3)S' = so that 
5' = or a = 3/2. 

If <S" = 0, the resulting solution is 

e' p = coe K/{1+K \ (3-41) 

e* = 30 (3-42) 
y/1 -w s$ 

S = s , (3-43) 

Mi = w sl M 2 = 0, (3-44) 

P 1 = KWi = Kw , (3-45) 

P 2 = W 2 = 0, (3-46) 

where cq, sq and wo are constants and 4K = wqSq{K 2 + 6K + 1). — 1 < if < 2y/2 — 3 
and < -RT < 1 must be satisfied for the positive energy density. With the coordinate 

transformations t' = c Sq 2K/{1+K) a~ 2K /^ 1+K ^ (1 - 2K/[a(l + K)])' 1 ^'^^) and 
r' = so r ) the solution is found to be a static solution: 

ds 2 = _ r ,4K/(l+K) dt ,2 + ^+ + 6 ^+ 1 d /2 + /2 ^ 2) (3 . 47) 
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4K 2 

8irGp = 8ttGKii = -— r— r. (3-49) 

F (K 2 + 6K + l)r' 2 v ' 

This solution is singular at the physical center r' = 0. We call this solution the 
singular static solution. It is noted that this solution is also of the first kindEj). 

For a = 3/2, equations (3-5) and (3-7) are degenerate. Equation (3-10) and 
3S'/S + [2/(1 + K)]W{/Wi = give 

Equations (3-50) and (3-17) give the evolution equation for S: 



2(1 + K)y' + (3 - K)y + 3(1 - K)(l + K)y 2 = 0, (3-51) 
where y = S'/S. This equation can be integrated as 



S = s 1 \l-s £ 2( 3 i+x) | 3(i-k) ( ( 3 . 52 ) 

where so and s\ are integration constants. An identity related to Mi can be obtained 
from equations (3-4), (3-6) and (3-8): 



WxS(S + S') + KW X SS' = 1 - e- 2q, {S + S') 2 . (3-53) 
Eliminating <£>' from equation (3T8) with using equations (3T7), we obtain 

2{S + S')(SS" + 2SS' + 5 /2 ) = 5'(1 + KiyiS^e 2 ^. (3-54) 
From equations (3-51), (3-53) and (|3-54j) , we obtain 

w _ (K 2 + 6K + 1)S 2 + (1 + K)(9K + 2)SS' + (3^ + !)(! + K) 2 S' 2 

e - (1 + K) 2 ' [6 ' bb) 

and 



w = Pi = K[A + 3(1 + K)S'/S] 

1 K (K 2 + QK + l)S 2 + (l + K)(9K + 2)SS' + (3K + l)(l + K) 2 S' 2 ' 

(3-56) 

exp(^), Mi and M2 can be obtained from equations (3-50), (3-8) and (3-9), respec- 
tively: 



12 



H. Maeda, T. Harada, H. Iguchi and N. Okuyama 



e* = c e K/{1+K) S 3K / 2 , (3-57) 
Ml = S[l — e~ 2>I '(S + S') 2 }, (3-58) 

M 2 = ^e~ 2 *SS' 2 , (3-59) 

where cq is an integration constant. Since 3S' /S+ [2/(l+K)] [W{/Wi] = gives oc 
W^ 2 ^ l+K \ this solution is singular at the physical center 5 = 0. For sq > 0, this 
solution represents collapsing shells, since S decreases as £ decreases from £ = +oo 
which means that i increases from t = +0 for constant r. A shell for some constant 
r collapses to a singularity 5 = at a finite time t\ = (2/3)s 3 ( 1+x ")/( 3 K \z/i p or 
so < 0, this solution represents expanding shells, since S increases as £ decreases 
from £ = +oo which means that t increases from t = +0 for constant r. In this case, 
S cannot be zero. We call this the dynamical solution (A). 

3.2. Self- similarity of the the zeroth kind 

In the case of self-similarity of the zeroth kind, the Einstein equations imply 
that the quantities fj,,p and m must be of the form 



2Gm = r[Mi(£) + r 2 M 2 (£)], (3-60) 
STrGfi = ^[Wi(0 + r 2 W 2 (0], (3-61) 

8irGp = ±[P 1 (0 + r 2 P 2 (0}, (3-62) 

where £ = re~*. A set of ordinary differential equations is obtained when one 
demands that the Einstein equations and the equations of motion for the matter 
field are satisfied for the O[r ] and 0[r 2 ] terms separately. The resulting equations 
for a perfect fluid (2-4)-(2-ll) reduce to the following: 



Mi + M[ 


= W^iS + S'), 


(3-63) 


3M 2 + M 2 


= W 2 S 2 (S + S'), 


(3-64) 


M[ 


= -P!S 2 S', 


(3-65) 


M' 2 


= -P 2 S 2 S', 


(3-66) 


Mi 


= S[l-e- 2 *(S + S') 2 ], 


(3-67) 


M 2 


= SS' 2 e- 2)p , 


(3-68) 


(Pi + W{)& 


= 2Pi - P{, 


(3-69) 


(P 2 + W 2 )& 


- P' 


(3-70) 


w{s 


= -(Pi + Wi)(<F , S + 25 / ), 


(3-71) 




= -(P 2 + W 2 )(<f , 1 S + 2 1 S / ), 


(3-72) 


s" + s' 


= SV + (5 + 5')!^, 


(3-73) 
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S'(S' + 2&'S) = W 2 S 2 e 2<p , (3-74) 

2S(S" + 2S') - 2&'S(S + S') = -S' 2 - S 2 + e 2 * (1 - W X S 2 ), (3-75) 

2S{S" - $'S') + S' 2 = -P 2 S 2 e 2<p , (3-76) 

(5 + S')(S + S' + 2&S) = (1 + PxS^e 2 *, (3-77) 



where the prime denotes the derivative with respect to ln^. In a vacuum case, the 
Minkowski spacetime can be obtained since Mi = M 2 = 0. 

3.2.1. EOS (I) and (II) 

Subtracting equation (3-77) from equation (3-75) and eliminating S" with using 
equation (3-73), we obtain 

2<2>' = (P 1 + WxY 2 * . (3-78) 
Then equations (3-69) and (3-70) result in 

e 2 ' I, (P 1 + W 1 ) 2 = AP 1 -2P{, (3-79) 
e 2 *{P 1 + W 1 ){P 2 + W 2 ) = -2P 2 '. (3-80) 

If a perfect fluid obeys EOS (I), we find from equations (3-61) and (3-62) that 

P 1 = W 1 = 0, P 2 = K{MG) l -^W^. [case (A)] (3-81) 
If a perfect fluid obeys EOS (II), we find from equations (3-61) and (3-62) that 

*- w >-°- ft - (8 ^-.^ ( y '-^i) T - |case(B)l (3 ' 82) 

Since Pi = W\ = 0, M\ = can be obtained from equations (3-63) and (3-65). 
Equation (3-67) gives (S + S') 2 = e 2 * . From equations (3-78) and (3-80), P l =W 1 = 
implies = and P 2 = Po, where po is a constant, which means that W 2 = wo, 
where wq is a constant. We set = cq, where Co is a constant. Then equation 
(3-72) gives (&S + 2S')(P 2 + W 2 ) = 0. P 2 + W 2 = means p + n = 0, which will 
be treated in the next subsection. If & f S + 2S' = 0, the equations (3-64), (3-66) and 
(3-68) give the evolution equation for S; 

;!(§-) -(Po + ^o)|r-wo = 0. (3-83) 

The solution is S = soC q where q is a constant, q = can be obtained from equations 
(S + S') 2 = e 2 ^ and V'S + 25' = 0, which means that S = s and P 2 = W 2 = 0. 
Therefore the spacetime must be vacuum. 
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3.2.2. EOS (III) 

If a perfect fluid obeys EOS (III), we find from equations (3-61) and (3-62) that 



Pi = KWi, P 2 = KW 2 . [case (C)] (3-84) 

W{W 2 = can be proved. Assuming that W\ 7^ and W 2 7^ 0, we find from 
equations (3-71) and (3-72) that W[/W x = W^/W 2 , while equations (3-69) and (3-70) 
give that 2 — W[/W\ = — W 2 /W 2 - These two equations contradict each other. Hence 
W\W 2 = can be concluded. 

In the case of W\ = P\ = and W 2 / 0, the discussion in the previous subsection 
applies and gives that 



2$ 2 

2 = Po, 



(3-85) 

e 2 * = {S + S'f, (3-86) 

Mi = 0, (3-87) 

M 2 = Cq 2 SS' 2 , (3-88) 

P 2 = -W 2 = po, (3-89) 

where cq and po are constants. The evolution equation for S is 

-Po, (3-90) 
with a general solution 



3 fS'^ 2 



S = 8 ^-M JL ^-. (3-91) 



This spacetime is the de-Sitter spacetime with the metric 



ds 2 = _ c 2^2 + e 2^^*[ (1 _^_ c 2 po/3) 2^ 

(3-92) 

With the coordinate transformations t' = c^t and r' = sor 1_ v /_c oPo/ 3 ) the solution 
in a more usual form 



^2 = _ dt !2 + e 2^%t' {dr ,2 + r ,2 dQ 2^ (3 . 93) 

2 Gm = -^A 3 ^', (3-94) 
8ttGp = -8ttG// = 8vrGp , (3-95) 



is obtained. 
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In the case of W2 = P2 = and W\ / 0, equations (3-64) and (3-66) give M2 = 0, 
and then S' = can be found from equation (3-68). If K = —1, Pi = — W\ = po 
is obtained from equation (3-71) where po is a constant. Equation (3-69) then gives 
Po = which means that the spacetime must be vacuum. For K 7^ — 1, the resulting 
solution is 



e * = CQ e K/{i+K \ 


(3-96) 


& so 
6 = 1 9' 


(3-97) 




S = s , 


(3-98) 


Mi = w sl M 2 = 0, 


(3-99) 


P 1 = KWi = Kw , 


(3-100) 


P 2 = W 2 = 0, 


(3-101) 



where Co, so and wo are constants and 4K = wqSq(K 2 +6K +1). With the coordinate 

transformations t' = cos Q 2K ^ 1+K \—(l + K) / (2K)]e~ * and r' = sor, this solution 
represents the same singular static solution as equations (3-41)-(3-46). — 1 < K < 
1\pl — 3 and < K < 1 must be satisfied for the positive energy density. 

3.3. Self- similarity of the infinite kind 

In the case of self-similarity of the infinite kind, the Einstein equations imply 
that the quantities [i,p and m must be of the form 



2Gm = Mi (OA 2 + M 2 (0, (3-102) 
8ttG/x = WiiO/t 2 + W 2 (0, (3-103) 
87rGp = P 1 (0/t 2 +P 2 (0, (3-104) 

where £ = r/t. A set of ordinary differential equations is obtained when one demands 
that the Einstein equations and the equations of motion for the matter field are 
satisfied for the O[t ] and Oft~ 2 ] terms separately. The resulting equations for a 
perfect fluid (2-4)-(2-12) reduce to the following: 



M[ 


= W X S 2 S', 


(3-105) 


M' 2 


= W 2 S 2 S'. 


(3-106) 


2Mi + M[ 


= -P!S 2 S', 


(3-107) 


M' 2 


= -P 2 S 2 S', 


(3-108) 


Mi 


= e~ 2 *SS' 2 , 


(3-109) 


M 2 


= S(l - e- 2 *S' 2 ), 


(3-110) 


(Pi + Wi)#' 


= ~P{, 


(3-Hl) 


(P 2 + w 2 )& 


= ~Pi 


(3-112) 
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(2Wi + W[)S 


= -(P 1 + W 1 )(V'S + 2S'), 


(3-113) 




= -(P 2 + W 2 )(^S + 2S'), 


(3-114) 


s" 


= S'(& + V), 


(3-115) 


W 1 S 2 e 2<p 


= S'(2W'S + S'), 


(3-116) 


(1 - W 2 S 2 )e 2 * 


= 2SS" + S' 2 - 2VS'S, 


(3-117) 


(l + P 2 S 2 )e 2 * 


= S'(2<P'S + S'), 


(3-118) 


-P 1 S 2 e 2 ^ 


= 2S"S + 2SS'{\-$') + S' 2 , 


(3-119) 


-P^e 2 * 


= S" + V"S + (1 - #' + V'){S' + &S), 


(3-120) 


-P 2 Se 2 * 


= -S" - <P"S + - &)(S' + &S), 


(3-121) 



where the prime denotes the derivative with respect to ln£. From equations (3-115), 
(3-116) and (3-119), 

(Pi + W^Se 2 * = -2S', (3-122) 
can be obtained, while from equations (3-115), (3-117) and (3-118), 

P 2 + W 2 = 0, (3-123) 

can be obtained. A vacuum case is not compatible since equation (3-122) gives 
S' = 0, which contradicts equation (3-118). 

3.3.1. EOS (I) and (II) 

If a perfect fluid obeys EOS (I), we find from equations (3-103) and (3-104) that 

P 1 = W 1 = 0, P 2 = K^Gf^W^. [case (A)] (3-124) 
If a perfect fluid obeys EOS (II), we find from equations (3-103) and (3-104) that 

Pl = H/l=0 ' P2= (faG^x (^-^l) 7 - [case (B)] (3-125) 

If Pi = W\ = 0, equation (3-123) implies that these are included in the case of 
p = Kji, which will be treated in the next subsection. 

3.3.2. EOS (III) 

If a perfect fluid obeys EOS (III), we find from equations (3-103) and (3-104) 

that 

P 1 = KW U P 2 = KW 2 . [case (C)] (3-126) 

In this case, equation (3-123) implies that K = —1 or P 2 = W 2 = 0. For K = — 1, 
equations result in 
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P 1 = Wi = 0, (3-127) 

P 2 = -W 2 = —j, (3-128) 

Mi = 0, (3-129) 
M 2 = S = s , (3-130) 

where sq is a constant. This constant satisfies 1 = K(8nG) 1 ~' y Sq 2 ^" 1 ^ for case (A), 
while —1 = if/[(87rG) 7-1 m2'][7/(7 — l)] 7 s 2 ^ 7 ^ for case (B). The metric functions 
^ and \P are obtained by the integration of the following differential equations: 

0" = _JL e 2 *-<£ /2 + #V, (3.131) 

4 

= - - 1)$'. (3-132) 

In this solution, the perfect fluid is a cosmological constant since p + fj, = 0. There is 
no physical center because of the constant circumferential radius. Although we have 
not obtained the analytic form of general solutions of equations (3-131) and (3-132), 
there is a special solution in which ^> is constant. For this case, the above equations 
are explicitly integrated as 



e # = A sin (^ ln ^) +£cos (^ ln ?) » (3-133) 
e* = Cl , (3-134) 



where A, B and c\ are constants. With coordinate transformations t Cl ^ s ° = t! and 
r ci/s _ r ' -j-^g so l u tion in another form 



ds 2 = -fo^o/ci-2 [ Agin (i nr '_l nt ') +jB cos (In/ -Int')] 2 dt' 2 +^{dr l2 +r' 2 dQ 2 ) 



-l ' 



(3-135) 



2Gm = s , (3-136) 
8ttGp = SirGfi = -\, (3-137) 

is obtained. We call this solution /1-cylinder solution. It is clear that this new 
solution is not identical with but closely related to the solution found by NariaiE$. 
For K / -1 and P 2 = W 2 = 0, equation (3-111), (3-113) and (3-115) give 



18 



H. Maeda, T. Harada, H. Iguchi and N. Okuyama 



It is noted that S' cannot be zero since S' = contradicts equation (3-118). Sub- 
stituting equations (3-116) and (3-119) for P\ = KW\, and eliminating and 
W[/W\ with using equations (3-111), (3-113) and (3-138), then we obtain the evolu- 
tion equation for S: 

2(1 + K)y' + 2(1 - K)y + 3(1 -K)(l + K)y 2 = 0, (3-139) 
where y = S'/S. From equations (3-139) and (3-116)-(3-119), we obtain 



e 



2& 



2K 
l + K 



SS' + (1 + 3K)S' 2 . (3-141) 



Wi, Mi and M 2 can be obtained from equation (3-138), (3-109) and (3-110), respec- 
tively: 

Wi = § = li/or 2 /* 1 -" ^- 1 ^- 2 , (3-142) 

Mi = e~ 2,p SS' 2 , (3-143) 

M 2 = m 2 = 5(1 - e~ 2 *S' 2 ), (3-144) 

where u>o and m 2 are constants. For K ^ 1, equation (3-139) can be integrated to 
give 



1 S , = si|l-s ?" T + 1? | 3(1 - if) , (3-145) 

where s$ and si are integration constants. Since equation (3-140) implies W± oc 
£-2/{i+K) s -2 S ,-i K _ So £_(i-jo/(i+A-)|-(i+A-)/(i-JO this solution is singular 

at the physical center S = 0. For so > 0, this solution represents collapsing shells, 
since S decreases as £ decreases from £ = +00 which means that i increases from 
i = +0 for constant r. A shell for some constant r collapses to a singularity S = 
at a finite time £1 = s ^"V. For so < 0, this solution represents expanding 

shells, since S increases as £ decreases from £ = +00 which means that t increases 
from t = +0 for constant r. In this case, S cannot be zero. We call this the dynamical 
solution (B). For K = 1, equation (3-139) can be integrated to give the flat FRW 
solution: 



e # = c , (3-146) 
^ = ^„r 1/S , (3-147) 
S = s C 1/3 , (3-148) 
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Mi = lc 2 s 3 C\ (3-149) 

M 2 = 0, (3-150) 

Pi = Wi = ^c - 2 , (3-151) 

P 2 = ^ 2 = 0, (3-152) 

where Co and so are constants. With the coordinate transformations t' = c$t and 
r ' = (W c c/ 3 ) r_1 ^ 3 i the solution in a more usual form 

ds 2 = -dt' 2 + t' 2 '\dr' 2 + r'W), (3-153) 

2Gm = (3-154) 

8ttG P = 8ttG^ = ^, (3-155) 



is obtained. 



§4. ^ parallel to 



4.1. Self- similarity of the second kind 
In this case, we can choose £ to be 



and then the metric can then be written as 



ds 2 = _ t 2(a-l) e 2<P(r)^2 + ^2 + (4 . 2) 

In this case, the Einstein equations imply that the quantities m, /j, and p must be of 
the form 



2Gm = tMi(r) + i 3 " 2Q M 2 (r), (4-3) 
8vrG^ = t~ 2 Wi{r) + ^ 2a VF 2 (r), (4-4) 
SnGp = r 2 Pi(r) + ^ 2o P 2 (r), (4-5) 

and a set of ordinary differential equations is obtained when one demands that the 
Einstein equations and the equations of motion for the matter field are satisfied for 
the O[t ] and 0[t 2 ~ 2a ] terms separately. The resulting equations for a perfect fluid 
(2-4)-(2-ll) reduce to the following: 



Mi = -PiS 3 , (4-6) 
(3-2a)M 2 = -P 2 S 3 , (4-7) 
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M[ 


= WiS 2 ^, 


(4-8) 


M' 2 


= W 2< S 2 S', 


(4-9) 


Mi 


= 5[l-5' 2 ], 


(4-10) 


M 2 


= S 3 e- 2 *, 


(4-11) 


3Pi 


= -Wi, 


(4-12) 


3P 2 


= (2a - 3)W 2 , 


(4-13) 


-Pi 


= (Pi + WiX 


(4-14) 


-PL 


= (P 2 + W 2 X 


(4-15) 





= 


(4-16) 


2S"S + S' 2 


= 1 - WiS 2 , 


(4-17) 


3 


= W 2 e 2 *, 


(4-18) 


S« + 2$>S'S 


= 1 + PiS 2 , 


(4-19) 


2a - 3 


= P2e 2 *, 


(4-20) 



where the prime denotes the derivative with respect to r. A vacuum spacetime is 
not compatible with this case since W 2 = contradicts equation (4-18). 

4.1.1. EOS (I) and (II) 

If a perfect fluid obeys EOS (I), we find from equations (4-4) and (4-5) that 

a = 7, P 1 = W 2 = 0, P 2 = ^JL— Wl, [case (A)] (4-21) 

or 

a = i P 2 = W 1= 0, Py= (87r g )7 _! ^2 7 - [case(B)] (4-22) 
If a perfect fluid obeys EOS (II), we find from equations (4-4) and (4-5) that 



a = 7, Pi = 0, P 2 = * = (7 - 1)W 2 , [case (C)] (4-23) 



(8irG)" / ~ 1 m2 

or 

1 K 



a = -,P 2 = 0, P 1= (87rG)7 -i m 7 ^2 7 = (7-l)^i- [case (D)] (4-24) 

Since W2 cannot be zero, case (A) is excluded. If Pi = 0, equation (4-12) gives 
W\ = 0, which implies W 2 = for case (C). Hence case (C) is also excluded. 

For P 2 = 0, equation (4-20) gives a = 3/2 which implies 7 = 2/3 for both cases 
(B) and (D), and then equations reduce to the following: 



e = c , 



(4-25) 
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'o, (4-26) 

(4-27) 

(4-28) 
(4-29) 
(4-30) 

where cq and po are constants. Pi = —W\/2> implies that the spacetime must be 
vacuum in case (B). We require that W\ is positive so that po < 0. The solution is 



Pi = 


-3^i 




3 


w 2 = 






C 


Mi = 


-PtS 3 


M 2 = 


Cq D , 


PoS 2 = 


1-S' 2 



s = 



sin y/—por, 



(4-31) 



which gives the metric 

ds 2 = -cltdt 2 + t 2 (dr 2 + — !— sin 2 yf^rdfl 2 ] . 

V (-Po) ' / 



(4-32) 



With the coordinate transformations t' = (2co/3)t 3//2 and r' = ^/— po r > this solution 
is found to be the closed FRW solution with dust and p = — fj,/3 comoving fluids: 



ds 2 
2Gm 
8nGp 



/ 1 \ 4 /3 

dt >2 + (_ po )-l ^_Lj t'4/3 + gin 2 r , dQ 2^ (4 . 33) 

sin 3 r', (4-34) 



-(^)^o" (-Po)- 1/2 t' 2 / 3 + c 2 (-po)- 3/2 



3^ 3 ^_ 2 / 3/ 

2, 



8vrG^ = -3p (^) t'- 4 / 3 + ^'- 2 . 



(4-35) 
(4-36) 



_3_^ - J/:! 
4.1.2. EOS (III) 

If a perfect fluid obeys EOS (III), we find from equations (4-4) and (4-5) that 



P l = KWi, P 2 = KW 2 . [case (E)] (4-37) 

In this case, equations (4-12) and (4-13) imply that a = 3/2 or W\ = 0. a = 3/2 
implies P 2 = W 2 = from equation (4-13) or (4-20). W 2 = contradicts equation 
(4-18). Hence W\ must be zero. The resulting solution is 



ds 2 = -t 2{a - l) cldt 2 + t 2 {dr 2 + r 2 dQ 2 ), 
P 1 = Wi = 0, 



(4-38) 
(4-39) 
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P 2 = Qa -ljW 2 = (2a - 3)co \ (440) 

Mi = 0, (4-41) 
M 2 = CqV, (4-42) 

where Co is a constant. This is the flat FRW solution. With the coordinate transfor- 
mations (co/a)t a = t' and (a/co) 1//Q r = r' , the solution in a more usual form 

ds 2 = -dt' 2 + t'i{dr' 2 + r /2 c^ 2 ), (4-43) 
2Gm = a^r' 3 , (4-44) 

8ttGp = 8vrG f-a - 1 ) n = 2a ~ 3 t'~ 2 , (4-45) 

is obtained. This solution exists for — 1 < K < 1 (K / 0,-1/3) which implies 
< a < 3 (a + 3/2,1). 

4.2. Self- similarity of the zeroth kind 

Equations (4-2)-(4-16) can be used for this case with a = 0. A vacuum spacetime 
is also not compatible with this case. 

4.2.1. EOS (I) and (II) 

If a perfect fluid obeys EOS (I), we find from equations (4-4) and (4-5) that 

P 1 = W 1 = 0, P 2 = K^Gf^W^. [case (A)] (4-46) 
If a perfect fluid obeys EOS (II), we find from equations (4-4) and (4-5) that 

Fi = w ''= ' ^^('''-ri)' [ca?e(B)1 (447) 

Equation (4-13) gives P 2 + W 2 = so that above two cases are included in the 
case of p = K/j,, which will be treated in the next subsection. 

4.2.2. EOS (III) 

If a perfect fluid obeys EOS (III), we find from equations (4-4) and (4-5) that 

P 1 = KW U P 2 = KW 2 . [case (C)] (4-48) 

In this case, equations (4-12) and (4-13) imply that W\ = P\ = since W 2 cannot 
be zero. The resulting solution is 



ds 2 = — -t~ 2 dt 2 + t 2 {dr 2 + r 2 dQ 2 ), (4-49) 

{-Po) 

P 1 = W 1 = 0, (4-50) 
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P 2 = -W 2 = po, (4-51) 
Mi = 0, (4-52) 

M 2 = ^V, (4-53) 

where po is a constant. This is the de-Sitter solution. With the coordinate transfor- 
mation t = exp(y/—po/3t'), the solution in a more usual form 



ds 2 = _ dt >2 + e 2^Jlt'^ r 2 + r 2 dQ 2^ (4 . 54) 

2Gm= ^r 3 e 3 ^', (4-55) 
8ttGp = SirGn = 8wGp , (4-56) 



is obtained. 

4.3. Self- similarity of the infinite kind 
In this case, we can choose £ to be 

?w - «£• (4 - 57) 

and then the metric can then be written as 

ds 2 = - e ™(r) dt 2 + dr 2 + S ( r ) 2 df2 2 . (4-58) 

In this case, the Einstein equations imply that the quantities m, \x and p must be of 
the form 



2Gm = M(r), (4-59) 
8vrG/i = W(r), (4-60) 
8vrGp = P(r), (4-61) 

and the Einstein equations and the equations of motion for the matter field (2-4)- 
(2-11) are written as 

M = S(l — S' 2 ), (4-62) 

M' = WS'S 2 , (4-63) 

(P + W)& = -P', (4-64) 

2SS" + S' 2 = 1 - WS 2 , (4-65) 

2<P'S'S + S' 2 = 1 + PS 2 , (4-66) 

where the prime denotes the derivative with respect to r. These equations give the 
Tolman-Oppenheimer-Volkoff equation. It implies that any spherically symmetric 
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static spacetime is a self-similar solution of the infinite kind in which kinematic 
self-similar vector is parallel to the fluid flow. In a vacuum case, the Schwarzschild 
solution can be obtained as 



M = m , (4-67) 
v = ±^/(S - m )S ± ^ In (s - ^ + m )S^j + c u (4-68) 

= co(l-^) 1/2 , (4-69) 



e 

where cq,c\ and mo are constants, while the Minkowski spacetime can be obtained 

as 

M = 0, S = r, e^ = c , (4-70) 

where Co is a constant. 

§5. ^ orthogonal to 17** 

5.1. Self- similarity of the second kind 
In this case, we can choose £ to be 

and then the metric can then be written as 

ds 2 = _ r 2a dt 2 + e 2^(t) dr 2 + ^VW. (5-2) 

In this case, the Einstein equations imply that the quantities m, \x and p must be of 
the form 

2Gm = rM^t) + r 3_2a M 2 (t), (5-3) 
8ttG// = r" 2 VFi(t) + r -2a W 2 (t), (5-4) 
8^Gp = r- 2 P x {t) + r- 2a P 2 (t), (5-5) 

and a set of ordinary differential equations is obtained when one demands that the 
Einstein equations and the equations of motion for the matter field are satisfied for 
the 0[r°] and 0[r 2 ~ 2a ] terms separately. It is noted that the solution is always 
singular at r = 0, which is correspond to the physical center. The equations for a 
perfect fluid (2-4)-(2-12) reduce to the following: 



Mi = WxS 3 , 



(5-6) 
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(3 - 2a)M 2 


= W 2 S 3 , 


(5-7) 


M[ 


= -P!S 2 S', 


(5-8) 


M' 2 


= -P 2 S 2 S', 


(5-9) 


Mi 


= S{l-e- 2 *S 2 ) 1 


(5-10) 


M 2 


= SS' 2 , 


(5-11) 


(2 - a)P 1 


= aWi, 


(5-12) 


aP 2 


= aW 2 , 


(5-13) 


W{S 


= -{P 1 + W 1 ){V'S + 2S'), 


(5-14) 


W 2 S 


= -{P 2 + W 2 ){V'S + 2S'), 


(5-15) 


(1 - a)S' 


= 


(5-16) 


-WxS 2 


= e -w S 2 - 1, 


(5-17) 


W 2 S 2 


= S' 2 + 2&'S'S, 


(5-18) 


-PiS 2 


= 1 - (1 + 2a)SV 2S \ 


(5-19) 


-P 2 S 2 


= 2S"S + S' 2 , 


(5-20) 


Pi 


2 -2* 

= a e , 


(5-21) 


-P 2 S 


= S" + W'S' + V"S + &' 2 S, 


(5-22) 



where the prime denotes the derivative with respect to t. A vacuum spacetime is 
not compatible with this case since Pi = contradicts equation (5-21). 

5.1.1. EOS (I) and (II) 

If a perfect fluid obeys EOS (I), we find from equations (5-4) and (5-5) that 

a = 7 , Pi = W 2 = 0, P 2 

or 

1 

a = -, P 2 = Wi = 0, Pi 

7 

If a perfect fluid obeys EOS (II), we find from equations (5-4) and (5-5) that 



(8vrG) 



h T W;<, [case (A)] (5-23) 



^§ypr W l [case(B)] (5-24) 



a = 7 , Pi = 0, P 2 = (87rG ^- lm 7 ^i 7 = (7 " 1)W 2 , [case (C)] (5-25) 

or 

a = i,P 2 = 0, Pi = (87rG) ^-i m 7 ^2 7 = (7-l)^i- [case (D)] (5-26) 

Cases (A) and (C) are excluded since Pi = contradicts equation (5-21). For 
P 2 = 0, equation (5-13) gives W 2 = 0, which results in Pi = for both cases (B) and 
(D), and therefore it can be concluded that there are no solutions in these cases. 
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5.1.2. EOS (III) 

If a perfect fluid obeys EOS (III), we find from equations (5-4) and (5-5) that 



P 1 = KW 1 , P 2 = KW 2 [case (E)]. (5-27) 

W1W2 = is concluded, since if not, equation (5-12) contradicts equation (5-13). 
Equation (5-21) implies that Pi cannot be zero, and hence W2 / so that W2 = 
P2 = can be concluded. Then equations (5-7) and (5-9) result in (3 — 2a)M2 = 
where M2 is constant. Equation (5-12) implies that K = a/(2 — a), and then 
— 1 < K < 1 implies that a cannot be 3/2. Hence M 2 = is concluded. The 
resulting solution is 



ds 2 


= -r 2a dt 2 


Mi 


= w sl, 


Pi 


a 

= - — 


2 — a 


P2 


= w 2 = 0, 



"—^dr 2 + sir 2 df2 2 , (5-28) 

1 - w s£ 

2 = 0, (5-29) 

= t^-^o, (5-30) 
I — a 

(5-31) 

where wq and sq are constants and (l + 2a — a 2 )vjQSQ = a(2— a). With the coordinate 
transformations Sq a t = t' and s^r = r', this solution can be written 



1 - W S% 

2Gm = w slr', (5-33) 

8,G P = SvrG-^-, = -^-^4. (5-34) 
2 — a 2 — a r z 

This solution represents the same singular static solution as equations (3-41)-(3-46) 
with K = a/ (2 — a). K = a/ (2 — a) gives K 7^ — 1 and the assumption a / 1 
corresponds to K 7^ 1. The positive energy density requires — 1 < K < 2\[2 — 3 and 
< K < 1 which correspond to a < (2\/2 - 3)/(\/2 - 1) and < a < 1. 

5.2. Self- similarity of the zeroth kind 

Equations (5-6)-(5-22) can be used for this case with a = 0. In a vacuum case, 
the Minkowski spacetime can be obtained since M\ = M2 = 0. 

5.2.1. EOS (I) and (II) 

If a perfect fluid obeys EOS (I), we find from equations (5-4) and (5-5) that 



P 1 = W 1 = 0, P 2 = K(87tG) 1_7 W 2 7 . [case (A)] (5-35) 
If a perfect fluid obeys EOS (II), we find from equations (5-4) and (5-5) that 
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p ' = w ' = ' p ' = ( 8jG * w fc-^r) 7 - [case(B)1 (5 ' 36) 

In this case, equations reduce to 



Mi 


= 0, 


(5-37) 


M 2 


= \w 2 s\ 


(5-38) 


e 2 ^ 


= s 2 , 


(5-39) 


VF 2 S 2 


= S' 2 + WS'S, 


(5-40) 


-P 2 S 2 


= 2S"S + S /2 , 


(5-41) 



which give the flat FRW metric 

ds 2 = -dt 2 + S 2 (dr 2 + r 2 dQ 2 ). (5-42) 
The evolution for S is governed by 

W 2 S 2 = 3S' 2 , (5-43) 
-P 2 S 2 = 2S"S + S' 2 . (5-44) 

These equations reduce to a evolution equation for S by using of EOS (I) or (II). S 
is generally not a power-law function of t in both cases. 

5.2.2. EOS (III) 

If a perfect fluid obeys EOS (III), we find from equations (5-4) and (5-5) that 

P 1 = KW U P 2 = KW 2 . [case (C)] (5-45) 
In this case, equation (5- 12) shows P\ = 0, and hence W\ = 0. Equations (5-43) and 

2 

(5-44) can be integrated to give S = sot 3(1+x) where so is a constant, which is the 

flat FRW spacetime for K / -1, while S = s e^ ( ~ Po)/3< (po = P* = -W 2 < 0), 
which is the de-Sitter spacetime for K = —1. With the coordinate transformation 
sqt = r', a more usual form of the flat FRW solution: 

ds 2 = -dt 2 + t^TK) ( dr « + r ^df2 2 ), (5-46) 
2Gm = (5-47) 

SirGp = SirGKfi = (5-48) 
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and the de-Sitter solution: 



2G m = tE^e^VF^V 3 , (5-50) 
3 

8-kGp = SirGn = p , (5-51) 



are obtained. 



5.3. Self- similarity of the infinite kind 
In this case, we can choose £ to be 

<"a? - 1- (5 ' 52) 

and then the metric can then be written as 

ds 2 = _ e 2r dt 2 + e 2Ht) dr 2 + ^(^2^2 (5.53) 

In this case, the Einstein equations imply that the quantities m, \x and p must be of 
the form 

2Gm = e~ 2r M x {t) + M 2 (t), (5-54) 
8ttG^ = e- 2r W"i(t) + W 2 (t), (5-55) 
8ttGp = e~ 2r Pi{t) + P 2 (t). (5-56) 

A set of ordinary differential equations is obtained when one demands that the 
Einstein equations and the equations of motion for the matter field are satisfied for 
the O[r ] and 0[e~ 2r ] terms separately. The resulting equations for a perfect fluid 
(2-4)-(2-ll) are 

S = M 2 = s , (5-57) 
M 1 = 0, (5-58) 
P 1 = W 1 = 0, (5-59) 

p 2 = -W 2 = -w = -\. (5-60) 

where sq and wq are constants. Equation (2-12) gives e 2 ^ = P 2 = —1/wq. Since 
e 2<? g-2 _ _y < 0, it can be concluded that there are no solutions in this case, 
independent of the form of the equation of state. 

§6. Summary and Discussions 



We have classified the kinematic self-similar perfect-fluid solutions with either 
EOS (I), (II) or (III). In most cases, the governing equations can be integrated to give 
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exact solutions, although there are a few exceptions. The analytic form of general 
solutions has not been obtained in the infinite-kind case with EOS (III) for K = — 1 
in which a kinematic self-similar vector is tilted. The results are summarized in 
table I-IV. We note that, independent of the form of the equation of state, any 
static solution is a kinematic self-similar solution of the infinite kind in the parallel 
case, since the governing equations give the Tolman-Oppenheimer-Volkov equation. 
It should also be noted that, independent of the form of the equation of state, 
kinematic self-similarity of the infinite kind in the orthogonal case is incompatible 
with a spherically symmetric spacetime. 

In the cases of EOS (I) and (II), i.e., the polytropic equation of state, the FRW 
solution is one of the compatible solutions. The closed FRW solution with dust 
and p = — /i/3 comoving fluids is the second- kind solution in the parallel case for 
EOS (II) with 7 = 2/3, while the flat FRW solution is the zeroth-kind solution in 
the orthogonal case for both EOS (I) and (II), in which the scale factor is not a 
power-law function of t in general. 

Next we summarize the case of EOS (III). The flat FRW solution is the second- 
kind solution in the tilted case for a ^ 3(1 + K)/2 and in the parallel case for 
a = 3(1 + K)/2, the zeroth-kind solution in the orthogonal case for —1 < K < 1 
and and the infinite-kind solution in the tilted case for K = 1. The de-Sitter 

solution is the zeroth-kind solution both in the tilted, parallel and orthogonal cases 
for K = — 1. The singular static solution, which is singular at the physical center, is 
both the second-kind solution for any a and zeroth-kind solution in the tilted case. 
Because of its staticity, this solution is also the infinite-kind solution in the parallel 
case. Positivity of energy density of the perfect fluid requires that -1 < K < 2^2-3 
or < K < 1 in these cases. The singular static solution is also the second-kind 
solution in the orthogonal case for K = a/ (2— a). Positivity of energy density implies 
that -1 < K < 2^2 - 3orO<K<lin this case. We have found two interesting 
exact solutions, which we call the dynamical solutions (A) and (B), respectively. 
The dynamical solution (A) is the second- kind solution in the tilted case for a = 3/2 
and K ^ —1, while the dynamical solution (B) is the infinite- kind solution in the 
tilted case for K 7^ —1,1. These solutions represent collapsing or expanding shells. 
We have also found another new exact solution, which we call /1-cylinder solution, 
yl-cylinder solution is the infinite- kind solution in the tilted case for K = — 1 and 
describes a spacetime with constant circumferential radius and no physical center. 
It should be noted that we have not obtained the analytic form of general solutions 
in the infinite-kind case for K = — 1 in which a kinematic self-similar vector is tilted. 
The dynamical solutions (A) and (B) and yl-cylinder solution could be important 
and we are now researching on these solutions ElP . 

In the vacuum case, the Schwarzschild solution is the second-kind solution in 
the tilted case for a = 3/2 and infinite-kind solution in the parallel case, while the 
Minkowski spacetime has kinematic self-similarity of the second kind in the tilted 
case for any a, of the zeroth kind in both the tilted and orthogonal cases and of the 
infinite kind in the parallel case. 

The polytropic perfect-fluid solutions compatible with kinematic self-similarity 
are the FRW solution and general static solutions. This result differs from that in 
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Newtonian case, for which there exists the polytropic Larson-Pens±on_aplution which 
is a self-similar solution describing a collapsing polytrope gasa^'t3''li3. The result 
in general relativity depends on the fact that Pi = or P2 = 0. The non-trivial 
solutions could arise when we assume the equation of state in which it is possible to 
have P1P2 7^ 0. We discuss the possible equations of state which are compatible with 
self-similarity in appendix |X[ A more general kind of self-similarity called partial 
homothety, which i mpo ses only equation (2-19) on a spacetime, could be compatible 
with EOS (I) or (II) & . Such cases will be investigated elsewhere. 

If the polytropic Larson-Penston solution is an attractor of generic non-self- 
similar collapse of polytropic gas in Newtonian gravity, as for the isothermal case, 
then how does the collapse proceed in the relativistic regime when we assume that 
the polytropic equation of state in relativistic regime is described by EOS (I) or 
(II)? For < 7 < 1, both EOS (I) and (II) are approximated by a dust fluid so that 
the generic collapse could converge to the spacetime whose central region can be 
described by the Tolman-Bondi solution^. For 1 < 7, EOS (II) is approximated by 
p = (7 — l)/i. For 1 < 7 < 1.036, there exists the general relativistic Larson-Penston 
solution, which is the self-similar solution of the first kind. Generic collapse, in 
the neighborhood of the center, could converge to this solution in an approach to a 
singularity. However, for 1.036 ^ 7, the attractor solution loses its attractive nature. 
For 4/3 < 7 < 2, the collapsing flat Friedmann solution, which is also a self-similar 
solution of the first kind, could be an attractor solution, since it has been shown 
by means of mode analyses that the solution has no unstable modes for spherical 
perturbations^. Stable self-similar solutions of the first kind for 1.036 ^$ 7 < 4/3 
have not been found so far. The candidate for an attractor solution is not known. 
For 2 < 7, the dominant energy condition could be violated in the central regime as 
the collapse proceeds, which is unphysical. 

In order to understand the whole picture of the generic collapse of a polytrope 
gas which obeys the equation of state (I) or (II) , full numerical simulations of grav- 
itational collapse will be quite helpful. 



Self-similarity 


Solution 


1st (tilted) 


none 


1st (||) 


none 


1st (_L) 


none 


2nd (tilted) 


none 


2nd (||) 


none 


2nd (_L) 


none 


zeroth (tilted) 


none 


zeroth (||) 


none 


zeroth (_L) 


flat FRW 


infinite (tilted) 


none 


infinite (| ) 


All static solutions with the EOS (1) 


infinite (_L) 


none 



Table I. Kinematic self-similar solutions for the EOS (I). || and _L denote the parallel and the 
orthogonal cases, respectively. See text for the values of parameters, a, K and 7. 
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Self-similarity 


Solution 


1st (tilted) 


none 


1st (||) 


none 


1st (_L) 


none 


2nd (tilted) 


none 


2nd (||) 


closed FRW 


2nd (_L) 


none 


zeroth (tilted) 


none 


zeroth (||) 


none 


zeroth (1.) 


flat FRW 


infinite (tilted) 


none 


infinite (||) 


All static solutions with the EOS (II) 


infinite (1.) 


none 



Table II. Kinematic self-similar solutions for the EOS (II). || and _L denote the parallel and the 
orthogonal cases, respectively. See text for the values of parameters, a, K and 7. 



Self-similarity Solution 



see Refs|2j),|l7|) 
FRW (see Ref [l0|)) 



Equation number 



1st (tilted) 



1st (||) 



1st (_L) 



see Ref 11) 



2nd (tilted) 



fiat FRW 
Singular static 
Dynamical solution (A) 



(3-41W3-46) 
(3-52), (pT55|)-(3-59) 



2nd (||) 



fiat FRW 



2nd (_L) 



Singular static 



(5-28)-(5-31) 



zeroth (tilted) 



de-Sitter 
Singular static 



(3-96)-(3-101) 



zeroth (|| 



de-Sitter 



zeroth (_L) 



fiat FRW 
de-Sitter 



infinite (tilted) 



flat FRW 

Dynamical solution (B) 
/1-cylinder solution 



(3-140H 3-145) 
(|3-135|)-(3-137) 



infinite (||) 
infinite (_L) 



All static solutions with the EOS (III) 



Table III. Kinematic self-similar solutions for the EOS (III). || and _L denote the parallel and the 
orthogonal cases, respectively. See text for the values of parameters, a, K. It is noted that the 
analytic form of general solutions has not been obtained in the infinite kind case with K = — 1 
in which a kinematic self-similar vector is tilted. 
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Self-similarity Solution 



2nd (tilted) 


Minkowski 




Schwarzschild 


2nd (||) 


none 


2nd (_L) 


none 


zeroth (tilted) 


Minkowski 


zeroth (||) 


none 


zeroth (±) 


Minkowski 


infinite (tilted) 


none 


infinite (||) 


Minkowski 




Schwarzschild 


infinite (_L) 


none 



Table IV. Kinematic self-similar solutions for the vacuum case. || and _L denote the parallel and 
the orthogonal cases, respectively. See text for the value of a parameter, a. 
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Appendix A 

Possible equations of state 

We have seen that the energy density and pressure are decomposed into one 
or two parts if the fluid is contained in a spherically symmetric spacetime which 
admit a kinematic self-similar vector. Here we consider what kind of equations of 
state is possible if the energy density and pressure are decomposed into the one 
or two parts. Since the following analysis based only on the form of the energy 
density and pressure, the equations of state may be restricted further if one demands 
that the matter content be a source of a spherically symmetric spacetime with the 
corresponding kinematic self-similar vector. 

We consider a barotropic equation of state, i.e., 

p = /(/*)■ (A-l) 

It implies that the pressure is a function only of the energy density. This class of 
equations of state is quite wide and it is often useful in many realistic and astro- 
physical situations. 

A.l. First kind 

For a tilted case, the energy density and pressure are expressed in the following 
form 

(A-2) 

P = (A-3) 

where £ = r/t. If we assume the barotropic equation of state, it is found that the 
function / must be of the form 

fix) = kx, (A-4) 
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where k is an arbitrary constant. This was proved by Cahill and Taub0). 

It can be easily proved that this result is also the case even for parallel and 
orthogonal cases. 

A. 2. Second kind 

For a tilted case, the energy density and pressure are decomposed into the fol- 
lowing form 

»=^ + f, (A-5) 
P-# + ^, <A-6) 
where £ = r/(at) 1 ^ a and a/0,1. The above can be rewritten as 

»-^ + ^ (A-7) 

where Wj = £~ 2 a~ 2 / Q VFi and Pi = £~ 2 a~ 2 ' a Pi. Hereafter we will omit the tildes 
for simplicity. 

Then, we can write down equation ( |A-1[ ) as 

t«/2 + i 2 " J I, f*/2 + f 2 J • lA '^ 

A partial differentiation with respect to t leads to 

a t(2/a)+l f3 J { t a/2 ^ t 2 )\ a t (2/a)+l t 3 J ' ^ W > 

where and hereafter in this section a prime denotes a derivative with respect to the 
argument. A partial differentiation with respect to £ leads to 

pm.m) _ ft {wm w 2 (o\ (wm wm\ , Ain 
i«72~ + -p - - 1 + — ; ■ + -12- j • (A-llj 

From the above two equations, we can eliminate /' as 

2 Pi(o n m)\ fwm + wm \ = ( 2 wm _ 2 wm \ ( m) + mo 



a t (2/ a )+i t 3 J \ t a / 2 t 2 J \ a t(2/«)+i i 3 J V i a / 2 t 2 

(A-12) 

Since a/0,1, the following relations must be satisfied. 

P 1 W 1 = W 1 P{, (A- 13) 

iW 2 = iy 2 p 2 , (A- 14) 

/W^ + aP 2 ^ = W X P' 2 + aW 2 .P[. (A-15) 
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For W±W 2 7^ 0, we can integrate equations (A- 13) and (A- 14) as 

Pi =k 1 W 1 , (A- 16) 

P 2 = k 2 W 2 , (A-17) 

respectively, where k\ and &2 are constants of integration. Substituting the above 
into equation (A- 15), we obtain 




(A-18) 



If k\ = k 2 = k is satisfied, we can find f(x) = kx. If not, we obtain 

W 2 = CW?, (A- 19) 

where C is a constant of integration. Then, we find 

kix + k 2 Cx a = f{x + Cx a ). (A-20) 

This equation of state is interpreted as a mixture of two comoving fluids with linear 
equations of state p\ = k±fj,i and p 2 = k 2 fi 2 with the relation fi 2 = Cfj,®, where 
p = p\ +p 2 and /U = Hi + fi 2 . This equation of state is also interpreted as a mixture 
of two comoving fluids with power-law equations of state p\ = kiC^ 1 ^ ^ a and 
p 2 = k 2 C a fj,2 with the relation /xi = Cjj^ ■ 

Next we consider the case in which I4 7 i=0. If W 2 7^ 0, then we obtain f(x) = 
kx + Cx l l a , where k and C are constants of integration. If W 2 = 0, Pi = P 2 = 
must be satisfied, which results in a vacuum spacetime. 

Finally we consider the case in which W 2 = 0. If W\ / 0, then we obtain 
f(x) = kx + Cx a , where k and C are constants of integration. 

The above results are summarized by a class of functions 

k lX + k 2 x a = f(dx + C 2 x a ), (A-21) 

where k\, k 2 , C\ and C 2 are arbitrary constants. It is seen that the above class of 
equations of state includes EOS (I) for k 2 = C\ = and k\C 2 7^ or for k\ = C 2 = 
and k 2 C x + 0, EOS (II) for k 2 = and k x C x C 2 / or for k x = and k 2 C x C 2 + 0, 
and EOS (III) for k x C 2 = k 2 d, (h,k 2 ) / (0,0) and (d,C 2 ) + (0,0). 

It can be easily proved that this result is also the case even for parallel and 
orthogonal cases. 

A. 3. Zeroth kind 

For a tilted case, the energy density and pressure are decomposed into the fol- 
lowing form 



(A-22) 
(A-23) 
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where £ = r /e* . The barotropic equation of state requires the following relations 

WiP[ = PxWi (A-24) 
W1P2 = P1W2. (A-25) 



If W\ 7^ 0, then we find 



Pi = kWi, (A-26) 
P 2 = fcW 2 + C, (A-27) 



where /c and C are constants of integration. Then, we obtain f{x) = kx + C. 

If W\ = 0, then Pi = or W 2 = 0. For W 2 = 0, W2 is constant, which implies 
Pi = and P 2 = const. For Pi = 0, we cannot specify the function / any further. 

The above results are summarized as follows: if C^p = C^/i = is satisfied, the 
equation of state cannot be specified any further from the two assumptions, while, 
if not, the equation of state is given by 

fix) = kx + C, (A-28) 

where k and C are arbitrary constants. 

It can be easily proved that this result is also the case even for parallel and 
orthogonal cases. 

A. 4. Infinite kind 

For a tilted case, the energy density and pressure are decomposed into the fol- 
lowing form 

V = ^^ + W 2 (0, (A-29) 
P=^ + p 2(0, (A-30) 

where £ = r/t. 

We can easily find that the result is completely the same as that for the zeroth 
kind. If C^p = C^n = is satisfied, the equation of state cannot be specified any 
further from the two assumptions, while, if not, the equation of state is given by 

fix) = kx + C, (A-31) 

where k and C are arbitrary constants. 

It can be easily proved that this result is also the case even for an orthogonal 
case. For a parallel case, since C^p = C^/J, = is satisfied, the equation of state 
cannot be specified any further. 
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